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When calls offered to a primary group of trunks find all of them 
busy, provisions are often made for these calls to overflow to other 
groups of trunks. Such traffic overflow systems have been of interest 
for a long time, but recently overflow systems that allow for some 
calls to be queued have been of importance. In this paper we analyze 
a traffic overflow system with queuing. The system consists of two 
groups, a primary and a secondary. We consider two cases which 
differ in the treatment of demands waiting in the primary queue. We 
adopt a novel analytical approach, which considerably reduces the 
dimensions of the problem and simplifies the calculation of various 
steady- state quantities of interest. Our theoretical results include 
expressions for the loss probabilities, the probability of overflow from 
the primary to the secondary, and the average waiting times in the 
queues. 

I. INTRODUCTION 

In this paper, a particular overflow system with queuing is analyzed. 
The system consists of two groups, a primary and a secondary, with n k 
servers and qk waiting spaces, which receive demands from indepen- 
dent Poisson sources S k with arrival rates X* > 0, k = 1 and 2, 
respectively, as depicted in Fig. 1. The service times of the demands 
are independent and exponentially distributed with mean service rate 
H > 0. If all n 2 servers in the secondary are busy when a demand from 
S 2 arrives, the demand is queued if one of the q 2 waiting spaces is 
available; otherwise, it is lost (blocked and cleared from the system). 
Demands waiting in the secondary queue enter service (in some 
prescribed order) as servers in the secondary become free. 

If all ni servers in the primary are busy when a demand from Si 
arrives and there is a free server in the secondary and no demands 
waiting in the secondary queue, the demand is served in the secondary. 
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Fig. 1 — Mean flow rates for an overflow system with queuing. 

If there are no free servers, then the demand is queued in the primary 
queue, if one of the q\ waiting spaces is available; otherwise, it is lost. 
Two different cases are considered for the treatment of demands 
waiting in the primary queue. In case I, a demand waiting in the 
primary queue may enter service either in the primary, when a server 
becomes free, or in the secondary, if a server becomes free and no 
demands are waiting in the secondary queue. In case II, no overflow is 
permitted from the primary queue, so that a demand in the primary 
queue must wait for a server in the primary to become free. 

The case q x = and q 2 = corresponds to no queuing. This problem 
has been analyzed by Brockmeyer 1 in the particular case X2 = 0, 
corresponding to no secondary source. The case A2 »* is a special case 
of two groups of servers with a common overflow group. The first 
group has n x servers and the second none, while the overflow group 
has 712 servers. This problem has been analyzed by Kosten, 2,3 but his 
analysis differs from ours, as is explained in the next section. 

Generally, let pa denote the steady-state probability that there are 
i demands in the primary and / demands in the secondary. These 
probabilities satisfy a set of generalized birth-and-death equations, 
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which take the form of partial difference equations connecting neigh- 
boring states. Rath 4 has solved these equations numerically in case I, 
by using a Gauss-Seidel iteration technique. Here we carry out an 
analysis that reduces the dimensions of the problem, which may be 
considerable in cases of interest. The basic technique is to separate 
variables in regions away from certain boundaries of the state space. 
This leads to various eigenvalue problems for the separation constant. 
The eigenvalues are roots of polynomial equations. The probabilities 
Pij are then represented in terms of the corresponding eigenfunctions. 
The constant coefficients in these representations are determined from 
the boundary conditions, and the normalization condition that the 
sum of the probabilities is unity. In general, these constants have to be 
determined numerically. In case II, additional states are accessible, 
since in this case demands wait in the primary queue even if a server 
is free in the secondary. In this additional region, the probabilities pij 
are expressed in terms of a fundamental solution of a partial difference 
equation. 

A different procedure for reducing the dimensions of the problem is 
that of Herzog, Woo, and Chandy. 13 Their procedure, if applied to the 
problem of this paper, would involve the recursive determination of a 
number of solutions of a partial difference equation, in order to express 
the probabilities pa in terms of their values on certain boundaries. 
Except in case II, we have to solve only ordinary difference equations 
in our approach, because of the separation of variables. In case II, as 
noted above, it is necessary to obtain only a single solution of a partial 
difference equation in a certain region. 

Various steady-state quantities which are of interest may be ex- 
pressed in terms of the probabilities py. The quantities include the 
loss (or blocking) probabilities, the probabilities that a demand is 
queued, the probability of overflow from the primary to the secondary, 
and the average waiting times in the queues. These quantities may be 
expressed directly in terms of the constant coefficients which occur in 
the representations for the probabilities p u . Thus the steady-state 
quantities of interest may be calculated directly, once the coefficients 
have been determined from the boundary and normalization condi- 
tions, without having to calculate the probabilities p» y . Here again the 
reduction in the dimensions of the problem is valuable. 

Only the theoretical results are presented in this paper. Numerical 
results will be presented in a forthcoming paper by Kaufman, Seery 
and Morrison. 5 The reduction in the dimensions of the problem by 
means of the analysis presented in this paper leads to a dense matrix, 
rather than the much larger sparse matrix for determining the proba- 
bilities py. Nevertheless, the approach presented here is amenable to 
numerical calculation for moderately large problems. The largest case 
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run so far using this approach has been for ni = 40 = n 2 and q\ = 10 
= q 2 , for case I. Kaufman 6 has used a numerical technique involving 
matrix separability (block diagonalization) for these problems and, in 
fact, is led to one of the same eigenvalue problems. In addition, 
Kaufman has obtained numerical solutions by means of successive 
over-relaxation techniques. For the example mentioned above, our 
results for the steady-state quantities of interest agree to 12 significant 
figures, although some of the probabilities ptj, which range over 30 
orders of magnitude, agree to fewer significant figures. 

An iterative solution of two-dimensional birth-and-death processes 
has been discussed by Brandwajn. 14 Although his numerical results for 
some specific problems indicate better execution times than with 
Gauss-Seidel iteration, the successive over-relaxation technique, with 
appropriate choice of parameter, would seem to be superior. 

The theoretical results are built up in stages, commencing with the 
analysis of the overflow problem without queuing in Section II. Next, 
the case in which there is only a secondary queue is investigated in 
Section III. This analysis pertains to both case I and case II, since the 
two cases differ only in the treatment of demands queued in the 
primary. The overflow problem in which there are both a primary and 
a secondary queue is investigated in Section IV for case I, and the 
steady-state quantities of interest are calculated in Section V. Another 
approach for case I, when there is only a primary queue, is presented 
in Section VI. The analysis of case II, when there is only a primary 
queue, is given in Section VII. The investigation of case II, when there 
is both a primary and a secondary queue, is commenced in Section 
VIII, and the steady-state quantities of interest are calculated in 
Section IX. 

Properties of the eigenfunctions which occur in the representations 
of the probabilities py t and which are related to Poisson-Charlier 7 ' 8 
and Chebyshev 9 polynomials, are given in Appendix A. Appendix B 
shows that the corresponding eigenvalues are distinct. Moreover, the 
eigenvalues, which are roots of polynomial equations, enjoy an inter- 
lacing property which facilitates their numerical evaluation. 

Another approach to the solution of the overflow problem, which is 
useful when q\ is large, or even infinite, as is the analysis in Section VI 
for case I with q 2 = 0, will be presented in a separate paper. Also, an 
alternate approach, which is useful when q 2 is large, or even infinite, 
is presented in another paper. 15 

II. OVERFLOW WITHOUT QUEUING 

We begin by considering an overflow system without queuing. The 
system consists of two groups, a primary and a secondary, with n* > 
1 servers, which receive demands from independent Poisson sources 
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Sk with arrival rates A* > 0, k = 1 and 2, respectively. The service 
times of the demands are independent and exponentially distributed 
with mean service rate [i > 0. If all n 2 servers in the secondary are busy 
when a demand from S 2 arrives, the demand is lost (blocked and 
cleared from the system). If all /ii servers in the primary are busy when 
a demand from S\ arrives, then the demand is served in the secondary, 
if one of the n 2 servers is free; otherwise, the demand is lost. 

The overflow system corresponding to the case of no secondary 
source (or A 2 = 0) has been analyzed by Brockmeyer. 1 Also, the system 
described above is a special case of two groups of servers with a 
common overflow group. The first group has ni servers and the second 
none, while the overflow group has n 2 servers. This problem has been 
analyzed by Kosten, 2,3 but his analysis differs from ours, as will be 
explained. 

Let pij denote the steady-state probability that there are i demands 
in the primary and j demands in the secondary. These probabilities 
satisfy a set of generalized birth-and-death equations, 10 which may be 
derived in a straightforward manner. We define the traffic intensities 

ai = Ai//x, a 2 = A2//1, (1) 

and let 8y denote the Kronecker delta; 6,7 = 1 if i = j, and 8y = if 
i ¥" j. Then it is found that 

[ai(l - 8 ini 8jn 2 ) + a 2 (l — 8j„ 2 ) + i +j]pij 

= ai(l — 8io)pi-ij+ (1 — 5/o)(ai8f n , + ozlpij-i 

+ (1 - 8 ini )(i + l)p i+ u + (1 - 8 Jn2 )U + DPij + i, (2) 

for < i < n\ , < j < n 2 . The normalization condition is 

i 1 Pu = 1. (3) 

,=0 y=o 

It is known from the theory of finite Markov chains 11 that there is 
a unique solution to the system (2) and (3), and p,> > 0. Once this 
solution has been obtained, the steady-state quantities of interest may 
be calculated. In particular, if Lk is the probability that a demand 
arriving from Sk is lost, then 

"i 

Li=p n „„ 2 , L 2 = X pi, n . 2 . (4) 

1=0 

Note that a demand from the primary source is lost only if all the 
servers in both the primary and the secondary are busy. Such a 
demand overflows to the secondary if all ni servers in the primary are 
busy, but there is a server free in the secondary, so that the probability 
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of overflow is 

/l2 = Y/>n liy . (5) 

7-0 

We now proceed to solve the system (2). Ifi^ni, the variables in 
(2) may be separated, and there are solutions of the form a,-/?,- where 

(ai + i + p)oti — ai(l — 5,o)a £ -i + (i + l)a,+i, (6) 

for < i < m — 1, and 

[a 2 (l - S Jn2 ) +j- pWj = «2d - 8jo)Pj-i + (1 - S jn2 )(j + l)/3, +1 , (7) 

for <j < n 2 , and p is a separation constant. The solution of (6) may 
be expressed in terms of Poisson-Charlier polynomials. 7,8 We here 
denote the solution of (6) for which cm = 1 by 8t(p, fli). The properties 
of s,-(p, a) which we will need are given in Appendix A. 

For <y < n 2 — 1, a solution of (7) is ft = «,-(— p, a 2 ), <y < n 2 . The 
boundary condition at / = n 2 imphes that 

(n 2 - p)s„ 2 (-p, a 2 ) = a 2 Sn 2 -i(-p, a 2 ). (8) 

With the help of (169) and (170) in Appendix A, this equation may be 
written in the form 

ps„ 2 (l-p, a 2 ) = 0. (9) 

But s^(l - p, a 2 ) is a polynomial in p of degree n 2 . It is shown in 
Appendix B that its zeros are positive and distinct, and we denote 
them by p m , m = 1, • • • , n 2 . Then, with p = 0, we have 

pmSr^il - pm, d 2 ) = 0, JTl = 0, • • • , U 2 . (10) 

It follows that we may represent the probabilities /»,> in the form 

"2 

Pa = 2 c m Si(p m , ai)sj(-p m , a 2 ), < i < tii, <y < n 2 , (11) 

m-0 

where the constants c m are to be determined. It remains to satisfy the 
boundary conditions corresponding to i = n\ in (2), as well as the 
normalization condition (3). Hence, for i = n\ and j ^ n 2 , we must 
have 

(ai + a 2 + ni +j)pn v j 

= aip ni -ij + (ai + a 2 ){l - 8 j0 )p ni j-i + (j + l)/V,+i, (12) 

for < j < n 2 — 1. If we substitute (11) into (12), it is found, after 
reduction with the help of the recurrence relations in Appendix A, that 
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"2 

2 c m [p m s ni (l + p m , ai)Sj(-p m , a 2 ) 

m=0 

+ ais„,(p m , ai)sj(-l - p m , a 2 )] = 0, (13) 

for0<y'</i 2 - 1. 

The boundary condition corresponding to i = /ii and j = n 2 in (2), 
namely, 

(/ii + n 2 )p ni ,n 2 = aip^-i.r^ + (ai + a^p^.^-i, (14) 

implies, after reduction and the use of (10), that 

£ c m [p m s„,(l + p m , ai)s^-i(l - p m , a 2 ) 

m-0 

+ aiS n ,(pm, a\)Sn^-\(-pm, d 2 )] - 0. (15) 

As is to be expected, this condition is redundant, as may be verified by 
summing (13) from /' = to n 2 — 1, and using (171) in Appendix A. 
Thus the coefficients c m are determined by (13) only to within a 
multiplicative constant, which is determined from the normalization 
condition (3). From (3), (10), (11), and (171), it follows that 

c s„,(l, ai)s^,(l, a 2 ) = 1. (16) 

In general, the constants c m have to be calculated numerically. In the 
limiting case a 2 = 0, for which p m = m, from (10) and (167), we were 
able to obtain an analytical expression for c m , and to reproduce the 
results of Brockmeyer, 1 but we omit the details. 

We now comment on the approach adopted by Kosten. 2,3 Whereas 
we obtained an eigenvalue problem by satisfying the boundary condi- 
tions at j = n 2 , but not at i = m, Kosten obtained an eigenvalue 
problem by satisfying the boundary conditions at i = n.\, but not at / 
= n 2 . He obtains the eigenvalue equation ys n| (l + y, at) = 0, with 
roots yo = and y r , r = 1, • • • , n\ . However, since the variables cannot 
be separated in (2) in the region < i < Hi, < j < n 2 — 1, the 
corresponding eigenfunctions, which occur in the representation of /?,>, 
contain two terms for r = 1, •■■, »i, rather than one as in (11). 
Moreover, the eigenfunction corresponding to yo = 0, which corre- 
sponds to the solution for n 2 = oo, is an infinite sum. Hence it would 
seem that Kosten's approach is preferable if n 2 is significantly larger 
than m, but that the representation (11) is more suitable if n 2 SS tii. 

Once the constants c m have been determined, the probabilities p,y 
may be calculated from (11). We note, in particular, that, with the help 
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of (10) and (16), and (167) and (171) in Appendix A, it follows that 
5 s,(0, ai) __ a[/i\ 

y_0 »n,Ui All y a y r \ 

r=0 

This is the well-known result for the steady-state probability that 
there are i demands in the primary. From (4), (5), and (11), with the 
help of (171), we obtain 

"2 

1*2 = £ C « S »l(l + P»»» «l) S n 2 ( _ P'n, «2) (18) 

771-0 



and 



"2 

/12 = £ c m s„j(p m , ai)sn 2 -i(l - p m , 02). (19) 

771-0 

Also, from (4), (5), and (17), we have 

r ^ n _ 8, 1 (0,ai) _ a?'/wi! 

Li + 7i2 = I Am./ " g n /, \ " "^ ( ' 

7=0 Sn,U. «U y fl 7y r | 

r=0 

We remark that the quantities L\ , L2 and /12 may be obtained directly, 
once the constants c m are known, without calculating the probabilities 



III. SECONDARY QUEUE ONLY 

We next consider the overflow problem in which there are 92 S 1 
waiting spaces in the secondary, but there is no primary queue. Now, 
if all n 2 servers in the secondary are busy when a demand from S2 
arrives, the demand is queued if one of the q 2 waiting spaces is 
available; otherwise, it is lost. Demands waiting in the secondary queue 
enter service (in some prescribed order) as servers in the secondary 
become free. If all ri\ servers in the primary are busy when a demand 
from S\ arrives, and if there is a free server in the secondary and no 
demands are waiting in the secondary queue, then the demand is 
served in the secondary. If there are no free servers, the demand is 
lost. 

Again, we let />,; denote the steady-state probability that there are 
i demands in the primary and j demands in the secondary, either in 
service or waiting. It is convenient to introduce the function 
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Then, with k 2 = ai 2 + q 2 , it is found that 

[ai(l - S^xy-zi,) + 02(1 - SjkJ + i + min (y, n 2 )]pij 

= ai(l - Sio)pi-u + (1 - S J o)(a 1 8in 1 Xn 2 -j + a 2 )pu-i 

+ (1 - Sinjd + l)pi+ij + (1 - 8 jk2 )rmn(j + 1, n 2 )p iJ+ i (22) 

for < 1 < Aii , < y < & 2 • The normalization condition is 

"i i Pij = 1. (23) 

,=0 y-0 

If i 7^ Mi, the variables in (22) may be separated, and there are 
solutions of the form a,/?/, where a, satisfies (6) for < 1 < Aii — 1 and 

[a 2 (l - Sjk.) + min(y, ai 2 ) - p]fij 

= a 2 (l - 8 j0 )Pj-i + (1 - 6> 2 )min(y + 1, /i 2 )ft + i (24) 

for < 7 < £2. As before, a, = Si(p, aO for < i < Aii. For < y < 
ai 2 — 1, (24) reduces to (7), and so there are solutions /?> proportional to 
Sj(—p, a 2 ) for <y < ai 2 . Also, for ai 2 <y < k 2 , 

[a 2 (l - fi,-* 2 ) + AI2 - p]0j = a 2 /8,-j + ai 2 (1 - 8jk 2 )fy + i. (25) 

The solution of (25) may be expressed in terms of Chebyshev polyno- 
mials of the second kind, 9 U,(x). It is convenient to define 

and 

<t>j(p) = Vkr-Ap) ~ **,-y-i<P)- (27) 

The properties of these functions which we will need are given in 
Appendix A. We note here, however, that Uo(x) = 1, U-i(x) = and 
<f>/, 2 (p) = 1. It follows from (25) and (177) that /? y is proportional to 
tjtjip) for Ai 2 - 1 <y : < k 2 . 
Consequently, we take 

n _ f Sj(-p, a 2 )«K(p)> 0<y<Ai 2 , , 9ft . 

" [ s„ 2 (-p, afefe(p), Ai 2 -l<y<£ 2 , U8 > 

where 

s n ,-i(-p, a>)4>n,(p) = Sn 2 (-p, a 2 )4hi 2 -i(p). (29) 

This equation may be written in the form 

p[s„ 2 (l - p, a>)*„ 2 (p) - s n ,-Al - P, a 2 )* 92 -,(p)] = 0. (30) 
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Note that this equation reduces to (9) if q 2 — 0. The expression in the 
square brackets in (30) is a polynomial in p of degree k 2 = n 2 + q 2 . 
Appendix B shows that its zeros are positive and distinct, and we 
denote them by p m , m = 1, • • • , k 2 . As before, we let po = 0. 

It follows that we may represent the probabilities py, for < i < /ii , 
in the form 

k 2 

2 C m Si(p m , ai)Sj(—p m , a 2 )<t>n 2 (prn), <j < U 2 , 

m=0 

Pi j={ (3D 

k 2 

£ CmSiipm, ai)Sn 2 (-p m , a 2 )<j>j(pm), /l 2 < > < Aj 2 , 

71=0 

where the coefficients c m are to be determined. Note that (31) reduces 
to (11) if q 2 = 0, since 4>* 2 (p) = 1. It remains to satisfy the boundary 
conditions corresponding to i = ri\ in (22), as well as the normalization 
condition (23). Hence, in addition to (12) for <y < n 2 — 1, we must 
have 

(a 2 + n\ + n 2 )pn v n 2 = aip^-i.n, + (a x + a 2 )p ni , n . 2 -i + n 2 p ni , n . 2 +i (32) 

and 

[a 2 (l - 8jk 2 ) + ni + n^p^j 

- aip ni -u + a 2 p ni j-i + n 2 (l - 8/*„)/>n,j+i (33) 

for n 2 + 1 <j<k 2 . 

If we substitute (31) into (12), we find, after reduction, that 

£ c m [p m s„,(l + p m , ai)sj{-p m , a 2 ) 

m=0 

+ ais ni (pm, ai)sj(-l - p m , a 2 )](j)n 2 (pm) = 0, (34) 

for <y < n 2 — 1. Also, from (33) we obtain 

k 2 

2 C m p m Sn,(l + Pm, ai)Sn 2 (-p m , a 2 )<f>j(p m ) = (35) 

m=0 

for zi 2 + 1 <j < k 2 , and from (32) we find that 

k 2 

X C m [p m S n ,(l + Pm, Cli)Sn 2 {—pm, Cl 2 ) 

m=0 

-ais ni (p m , ai)Sn 2 -i{-p m , a 2 )]<$>n 2 {pm) = 0. (36) 
This last condition is redundant, as may be verified by summing (34) 
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from j = to n.2 — 1 and (35) fromy = n 2 + 1 to k 2 , and using the fact 
that p m satisfies (30). Hence, the coefficients c m are determined by (34) 
and (35) only to within a multiplicative constant, which is determined 
from the normalization condition (23). From (23) and (31), with the 
help of (27), (30), and (170) and (171) of Appendix A, it is found that 

cos„,(l, aOts^U, a 2 )* 92 (0) - s„ 2 -i(l, a 2 )^ 2 -i(0)] = 1. (37) 
We note, from (167), (168), and (176), that 

Ml, a) - 2 ^, ¥,(0) =i(-). (38) 

Once the constants c m have been determined, the probabilities py 
may be calculated from (31). In particular, with the help of (27), (30), 
(37), (167), and (170), it is found that 

. £ s,(0, oi) qj/t! 

r=0 

as is to be expected. Moreover, the steady-state quantities of interest 
may be calculated. If Lk is the probability that a demand from Sk is 
lost, then 

la- 2 />».;> L 2= 2 p.-, * 2 . (40) 

j—n 2 1=0 

Also, the probability of overflow /12 is as given by (5), and so, from 
(39), 

Li + 7 12 = 2 ft,j = Q > —5 • (41) 

S „,(l, a,) 2 oJ/H 

r=0 

From (31), it is found that 
*s 
L 2 = 2 c "> s «,(l + Pm, ai)s n .,(-pm, a 2 ), (42) 

m— 

since <£*.,( p) — 1. from (27), and 

*2 

7l2 = 2 C mS n ,(pm, ai)S„,-\(l - p m , 2 ) <j>n 2 (pm)- (43) 

m=0 

Other quantities of interest are the probability that a demand from 
the secondary source is queued on arrival, namely, 

Q 2 =2 IPv, (44) 

i=0 j=n ., 
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and the average number of demands in the secondary queue, 

V 2 =X £ U-n 2 )pij. (45) 

j=oy=n 2 +i 

But, if we let 

rj-lpij, (46) 

and sum on i in (22), we obtain 

[a 2 (l - 8jk 2 ) + n 2 ]rj = a 2 r,-i + n 2 (l - 8jk 2 )rj+i, (47) 

for n 2 + 1 <y < A 2 . It follows that 

n 2 rj = azTj-i, n 2 + 1 <J < £ 2 - (48) 

Hence, since L 2 = r* 2 , from (40) and (46), 

0=(— ) 2 ^2, n 2 ^j<k 2 . (49) 

From (45), (46), and (49), we obtain 

f-A.fe) 1 * <50) 



where 



Similarly, from (44), 



(51) 



Q 2 = T 92 RJ L 2 , (52) 

where 

r,(0 = V r* = 9 + (* - 1)4,(1). (53) 

The departure rate from the secondary queue is 

n t k 2 k 2 

R 22 = n 2 fi'Z 2 Pv = "2/i 2 r y> (54) 

i=0 j=n 2 +\ j=n 2 +l 

from (46). It may be verified from (49) that R 22 = X 2 Q 2 . This result 
holds since, in the steady state, the departure rate from the secondary 
queue is equal to the arrival rate to it. 

Now, according to Little's theorem, 10 the average number of de- 
mands in a queuing system is equal to the average rate of arrival of 
demands to that system times the average time spent in that system. 
If we apply this result to the secondary queue, we find that the average 
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waiting time of the demands which are queued in the secondary is 

W 2 =r-^-, (55) 

A2V2 

independently of the order of service within the queue. Now the 
average number of demands being served by the secondary group is 

X 2 = S t min(j, n 2 )pij. (56) 

,=0 y=0 

Hence, if we apply Little's theorem to the secondary group of servers, 
we obtain 

a 2 (l - L 2 ) + aj i2 = X 2 , (57) 

from (1), since the mean service rate is ju. From (31) and (56), with the 
help of (23) and (171) and (172) of Appendix A, it is found that 

X 2 = n 2 - £ c m s ni (l + p m , ai)sn 2 -i(2 - p m , a 2 )<j>n 2 (p m ). (58) 

The relationship (57) may be verified from (42), (43), and (58), with 
the help of the boundary conditions (34) and (35), although the 
verification is difficult. Consequently, the relationship (57) provides a 
useful check on the numerical calculation of the constants c m . 

IV. PRIMARY AND SECONDARY QUEUES: CASE I 

We now begin the investigation of the overflow problem in which 
there is a primary queue, with q\ > 1 waiting spaces, as well as a 
secondary queue, with q 2 waiting spaces. As before, if all n 2 servers in 
the secondary are busy when a demand from S 2 arrives, the demand is 
queued if one of the q 2 waiting spaces is available; otherwise, it is lost. 
Demands waiting in the secondary queue enter service as servers in 
the secondary become free. If all ri\ servers in the primary are busy 
when a demand from Si arrives, and if there is a free server in the 
secondary and no demands are waiting in the secondary queue, then 
the demand is served in the secondary. If there are no free servers, 
then the demand is queued in the primary queue, if one of the q\ 
waiting spaces is available; otherwise, it is lost. 

We will consider two different cases for the treatment of demands 
that are waiting in the primary queue. In the first case, to be considered 
now, a demand waiting in the primary queue may enter service either 
in the primary, when a server becomes free, or in the secondary, if a 
server becomes free and there are no demands waiting in the secondary 
queue. Thus overflow from the primary queue is permitted in the first 
case. In the second case, to be considered later, no overflow is permitted 
from the primary queue, so that a demand in the primary queue must 
wait for a server in the primary to become free. 

OVERFLOW PROBLEMS WITH QUEUING 1439 



We let 

ki = ni + q i} k 2 = n 2 + q 2 . (59) 

Then, in the first case, it is found that the steady-state probabilities/?,, 
satisfy the equations 

[ai(l - S iki Xj-n 2 ) + a 2 (l - Sjk 2 ) + min(i, ni) -I- min(>, n 2 )]pij 

= (1 - Xi-i-ntXii-i-Jiad 1 ~ 8 io)Pi-hj 

+ (1 - 8jk 2 )mhiU + 1, n 2 )pi, j+ i\ 

+ (1 - 5 /0 )[ai5,„ 1 Xn 2 -> + a 2 (l - X'-i-»iX«2-./)]Pv-i 

+ (1 - 5,a,)[(1 - Xi- ni Xn2-w)min(t + 1, n,) + n 2 Xi-n&n 2 ]pi+ij, 

(60) 

for < i < Ai, <y < k 2 . These equations were constructed so as to 
imply that 

Pu = 0, nx + l<i<k u </' < n 2 - 1, (61) 

since it is impossible for demands to be waiting in the primary queue 
when there is a server free in the secondary. The normalization 
condition is 

I 1PU+ I I A7=l- (62) 

i'=0 7=0 i'=n, + l >=n 2 

We assume for the time being that q 2 > 1. For < i < n } - 1 and 
</' < k 2 , (60) is equivalent to (22), and hence there are solutions of the 
form eufij, where a, = «(p, ai) for < i S m, j8> is given by (28), and p 
is a root of equation (30). For m + 1 < i < £i and n 2 + 1 <j < & 2 , the 
variables in (60) may be separated, and there are solutions of the form 
ytSj, where 

[aid - «.-*,) + n x + p]y« 

= ai y,_i + ni(l - 5,a,)y,+i, /ii + 1 < i < £i (63) 

and 

[a 2 (l - 8 jk2 ) + n 2 - p]8j 

= a 2 8j-i + n 2 (l - 8jk 2 )8 j+1 , n 2 + 1 <j < k 2 , (64) 

and v is a separation constant. Corresponding to (26) and (27), we 
define 

QM - feY" tf,(*+£+i) (66) 

W \ 2va 1 n 1 / 

and 

ft-(i>) = Qi.-itf) - fVt-i('). (66) 
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Then y, is proportional to 0,-(i») for ni < i < k\ and 5, is proportional to 
4>j(p) for n 2 <j<k 2 . We note that 0*, (?) = 1. 

Now 0„, (?) is a polynomial of degree q\ in ?. It is shown in Appendix 
B that its zeros are negative and distinct, and we denote them by *v, 
«?= 1, ...,£i. Then 

M"')=°> ^=1, •••,91- (67) 

Since the roots p m , m = 0, • • • , £ 2 of (30) are nonnegative and distinct, 
we have a joint set of distinct roots. In view of the above results, we 
may represent the probabilities pij in the form 

k 2 

£ CmSiipm, di)Sj(-p m , a 2 )6 n] (Pm)(}>n 2 (pm), </' < n 2 , 

m=0 

Pij= (68) 

* 2 

2 C m Si(p m , di)s„A-p m , a 2 )9n x (pmHj(pm), R 2 < j <: k 2 , 
m=0 

for < i < ni , and 

Py™ £ c m s„,(p m , ai)s„ 2 (-pm, a 2 )di(p m )<f>j(p m ) 

m=0 

+ 2 bAMhM, m < i < Ai, /i2 <y < fe, (69) 

<r=l 

where the constants c m and £v are to be determined. Note that the 
representations in (68) and (69) agree for i = n\, in view of (67). We 
also note that (68) reduces to (31) if q\ = 0, since 0*, (?) ■ 1. 

It remains to satisfy the boundary conditions at i = ni , < y < k 2 
and aty = n 2 , ni + 1 < i < Ai, as well as the normalization condition 
(62). Hence, in addition to (12) for <y < n 2 - 1, we must have 

(ai + a 2 + ni + n 2 )p ni ,n 2 

= OiPm-1^2 + (fll + a2)p»,^ 2 -l + («1 + «2)Pn, + l,n 2 + tt2Pn,,n 2 + l (70) 

and 

[ai + a 2 (l - 6> 2 ) + ni + n 2 ]p ni j 

= aiPm-ij + a 2 p ni j-i + raip« 1+ i,y + n 2 (l - Sjk 2 )p ni j+i (71) 
for n 2 + 1 Sj ^ A 2 . Also, 
[ai(l - 5,a,) + a 2 + ni + /i 2 ]p,.„ 2 

= aip,-i,„ 2 + («i + n 2 )(l - 5i*i)Pi+i,n 2 + n 2 pi,n 2 +\ (72) 
for Hi + 1 < i < Ai . 
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If we substitute (68) into (12), we find, after reduction, that 

*2 

£ c m [p m s„,(l + p m , ai)sj(-p m , a 2 ) 

+ aiS n] (pm, ai)Sj(- 1 - p m , a2)]0 ni (Pm)^n 2 (Pm) = 0, (73) 

for s; < » 2 - 1. Next, from (68), (69), and (71), with the help of 
(184), we obtain 

*2 

2 C m p m Sn 2 (—pm, a2)[S„,(l + p m , ai)0n,(Pm) 
m=0 

+ S ni (p w , ai)fi 9| -l(pm)]<f»/(pm) = »1 2 Mn,+l(iv)$/(iv), (74) 

<r-i 

for «2 + 1 ^y ^ ^2. Also, from (72), with the help of the recurrence 
relation (63) satisfied by 0i(v), it is found that 

*2 

Yt CmSn^pm, Oi)S„ 2 (— p m , O.2)[a20i(pm)<t>n 2 -l(Pm) 

m=0 

— /l 2 (l - &A,)0i+l(p;n)<J>n 2 (pm)] 

+ S Wi("^ r i^) " "2(1 - 6\a 1 )0, + i(jv)<M»v)] = 0» (75) 

for ni + 1 < i < k\ . It may be verified that the boundary condition (70) 
is redundant. 

The coefficients c m and bt are determined by (73) to (75) only to 
within a multiplicative constant, which is determined from the nor- 
malization condition (62). But, from (68), with the help of (27), (30), 
and (170) and (171) of Appendix A, it is found that 

*2 

£ Pij = CoSi(0, ai)0 ni (O)|>„ 2 (l, a 2 )* 92 (0) - s„ 2 -i(l, a 2 )¥ 92 -.(0)] (76) 

y'-o 

for < i < in. Then, from (62) and (69), with the help of (66) and 
(171), it follows that 

c s„,(l, ai)0 ni (O)[s„ 2 (l, a 2 )% 2 (0) - s„ 2 -i(l, a 2 )^ 2 -i(0)] 

*2 

+ Yi CmSn y (Pm, Ol)S„ 2 (— p m , Ct.2)Slq x -l{pm)tyq 2 (pm) 
m—0 

+ I .MVl<*>**<*> = 1. 

t-\ (77) 

Note that (77) reduces to (37) if q x = 0, since Q_i(i») = and 8 ki (v) = 1. 
Once the constants c m and b e have been determined, the steady-state 
probabilities p;y may be calculated from (68) and (69). We remark that 
the number of constants to be determined is only fe + 1 + q u whereas 
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the number of probabilities p„ is (ni + l)(k 2 + 1) + qi(q2 + 1), which 
in general is considerably larger. 

V. SOME STEADY-STATE QUANTITIES: CASE I 

Other steady-state quantities of interest are depicted in Fig. 1, which 
indicates the mean flow rates. The loss probabilities L x and L 2 are 
given by 

U = 1 Pk vh U=i p ilk2 , (78) 

/-fig i-o 

and the probabilities that a demand from the primary, or secondary, 
source is queued on arrival are 

Qi = £ I PU, Q2 - S .2 Ay- (79) 

<— »i /— nj '=0 y— Bj 

The probability that a demand arriving from the primary source 
overflows immediately is Ji 2 , as given by (5). 

Since the mean service rate is n, the mean departure rate from the 
primary queue to the primary servers is 

Ru - n lt i J J p* (80) 

i— «i+i y=«2 

while the mean rate of overflow from the primary queue to the 
secondary servers is 

*i 
Ru = n 2 fi £ pt,^. (81) 

The mean departure rate from the secondary queue is 

*1 *2 

R22 = n 2 n £ £ P'J- ( 82 ) 

It may be verified from (60) that 

fin + R12 - A1Q1, R 2 2 = \2Q2. (83) 

The average number of demands in the primary and secondary 
queues are 

Vi- 2 E (J-»ilW, V 2 -E S U-n2)pu. (84) 

t=n,+l >=/i2 «=0 y'=/«2 + l 

Also, the average number of demands in service in the two groups are 
Xi-2 E** + »H 2 2 P.; < 85 > 

i-O y=0 i=n, + l j=n. z 
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and 



71] 7l2 _ 1 k\ k 2 

,=0 >=0 i=0 j=nz 



(86) 



If we apply Little's theorem to the primary and secondary queues, we 
find that the average waiting times of the demands which are queued 
in the primary or in the secondary are given by 



Vi V 2 

AiQi A2Q2 



(87) 



respectively, independently of the order of service within each queue. 
Also, if we apply Little's theorem to the primary and secondary groups 
of servers, we obtain 

\i(l -In- In) -Rn-tiXi, \ 2 (1 ~ L 2 ) + \Ji 2 + Ri 2 = (jlX 2 , (88) 

since the mean service rate is ju. 

The steady-state quantities of interest may be expressed in terms of 
the constants c m and b ( by use of the representations in (68) and (69). 
From (78), with the help of (27), since 6 kl (v) = 1, it is found that 

Li = £ c m s ni (p m , ai)s„ 2 (-pm, a 2 )^t q2 (p m ) + £ &/*«*(*/)> (89) 



m=0 



/■=! 



and, with the help of (66) and (171), since <f>* 2 (p) ■ 1, it is found that 

L 2 = Yt CmSnzi-pm, a 2 )[s ni (l + pm, ai)8 ni (pm) 
m-0 

+ Sn 1 (p m ,a 1 )^ Ql -i(p m )]+ i bM q ,-i(ve). (90) 
Similarly, it follows from (79) that 

*2 

Ql = £ C m S ni (Pm, a^S^-pm, a 2 )[to gj (p m ) - l]¥ 92 (p m ) 



+ 2 ^(r/) - l]*ft(*). (91) 
Next, from (5), with the help of (171), it is found that 

*2 

/l2 = 2 C m S ni (Pm, ai)Sn 2 -l(l ~ Pm, a 2 )6 ni {pm)<t>n 2 (pm). (92) 

ro=0 

Also, from (80) and (81), it follows that 



R\i = fiifi 



"2 

2 C m S nj (pm, ai)Sn 2 (-pm, a 2 )Q Ql -l(pm)¥ V2 (pm) 



m=0 



+ J 6AriW*ft('/) 



/=i 



, (93) 
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and 

R12 = n 2 n 



£ C m S ni (p m , ai)Sn 2 (—pm, a 2 )Sl qi -i(pm)<f>n 2 (pm) 

71=0 

91 
+ X bsti gi -l{l>t)<t>n 2 (t>t) 



. (94) 



The first relationship in (83) provides a useful numerical check, in view 
of the expressions in (91), (93), and (94). 
From (84), with the help of (27), (182), and (183), it is found that 

— 

* 2 Cm 
— Oi J] — S„,(p m , ai)Sn 2 (-p m , a 2 ) , I r 92 (pm) 
m = l Pm 

ni 



1 -fl 9l (p m ) +— Qu-i(Pm) 



-a, 2 -¥«(*) 

•=1 vt 



Hi 



l-0gj(i'/)+ — 0,,-l(F/) 



(95) 



where A 9 (£) is as defined in (51). From (85), with the help of (62), (76), 
and (172), it follows that 

Xi - n\ - c s ni -i(2, ai)0„,(O) 

.[ft*(l, a 2 )^ 92 (0) - s„ 2 _,(l, <fe)¥ Vl -i(0)]. (96) 
Similarly, from (86), with the help of (62), (68), (171), and (172), 
X 2 = n 2 - X c m s„,(l + p m , ai)s„ 2 -i(2 — p m , a 2 )^ n ,(pm)<|>n 2 (pm). (97) 

m-0 

In view of (89), (90), (92), (94), (96), and (97), the relationships in (88) 
provide a useful numerical check. 
If we now define 



- I P«» 



(98) 



and sum on i in (60), for n 2 + 1 </' < A 2 , we obtain (47), and hence 
(48). Also, L 2 = r* 2 , from (78) and (98), and it follows that (49) holds 
and, from (79) and (84), that (50) and (52) hold. Moreover, from (82) 
and (98), 

R 22 = n 2 n £ r„ (99) 

/-JI2+] 

and it may be verified, as before, that R 22 = \ 2 Q 2 . This completes the 
calculation of expressions for the steady-state quantities of interest. 
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We have assumed that q 2 ^ 1. We remark, however, that the 
representations in (68) and (69) remain valid in the case q 2 = 0, so that 
k 2 = /i2 and foGO = 1. We have verified that (73) and (75) to (77) hold 
for q 2 = 0, and note that * (p) = 1 and ¥-i(p) = 0. We have also 
verified that (89) to (97) hold for q 2 = 0. Moreover, V 2 = and Q 2 = 0, 
which is consistent with (50) and (52), as is seen from (51) and (53). In 
the next section, we investigate the case q 2 = in a different manner. 

VI. ANOTHER APPROACH: CASE I 

We now give special consideration to the case q\ > 1, g 2 = 0, so that 
k 2 = n 2 . Then (60) for < i < »i - 1 and < j < n 2 leads to the 
representation 

Pij = 1 d m Si(p m , ai)sA-p m , a 2 ), 0< i< »,, 0<y< n 2 , (100) 

m=0 

corresponding to (11), where p m , m = 0, • • • , n 2 are the roots of (10), 
with po = 0. The coefficients d m are to be determined. The only other 
nonzero steady-state probabilities are/>,,„ 2 for /ii + 1 < i < Ai. But, 
from (60), with k 2 = n 2 , 

[ai(l - 5,-a,) + ni + /i2]p«,n 2 

= aipi-i,n 2 + (ni + n 2 )(l - fi,-*,)pi+i.« 2 ( 101 ) 
for ni + 1 < i < ^! . It follows that 

(ni + n 2 )pi,n 2 = aipi-i,n 2 , ni + l<i<ki. (102) 

Hence, 

A> 2 =(— T— ) Pn lin2 , Til <*<*!, (103) 

\7li + /l2/ 

and pn^ is given by (100). We remark that (103) holds if qi = oo, 
provided that the stability condition ai < /ii + n 2 is satisfied. 

It remains to satisfy the boundary conditions at i = in, <y < n 2 , 
as well as the normalization condition (62) with A 2 = n 2 . But (60) 
imphes that (12) must hold for <y < n 2 - 1, and from (100) it follows 
that, corresponding to (13), 



"2 



2 rfm[PmS ni (l + pm, Cli)Sj(-p m , Gt 2 ) 



m— 



+ ais ni (p m , ai)sj(-l - p m , a 2 )] = 0, (104) 

for < j < n 2 - 1. It may be verified that the boundary condition 
corresponding to i = ni and/ = n 2 , namely, 

= fli/>„,-i,^ + (ai + a 2 )/)„„ n2 -i + (ni + Ai2)jD ni +i,n 2> (105) 
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is redundant. In fact, (105) reduces to (14), which was shown to be 
redundant, since 

(/ii + n 2 )pn 1 +i.n 2 = aiPm.^ (106) 

from (103). 

The coefficients d m are determined by (104) only to within a multi- 
plicative constant, which is determined by the normalization condition. 
From (62), with k 2 = n 2 , it is found, with the help of (10), (53), (100), 
(103), and (171), that 

^a^o..> + (^"\(s±s) 

n 2 

• £ d m S ni (pm, d^S^-pm, d 2 ) = 1. (107) 

m-0 

We note that, for a\ < n\ + n 2 , 

, im (^-Y'\(hi±*l) _ — « — < 108) 

9 ,-oo \rii + n 2 ) \ a x ) (/ii + n 2 -ai) 

Once the coefficients d m have been determined from (104) and (107), 
the steady-state quantities of interest may be calculated with the help 
of (100) and (103). Since the calculation of these quantities is a 
straightforward extension of previous calculations, in view of the 
simplicity of the result in (103), we do not present the results here. We 
remark that the approach adopted here in the case q 2 = requires the 
determination of only n 2 + 1 constants d m , m = 0, • • • , n 2 , whereas the 
approach adopted in the two previous sections requires the determi- 
nation of the n 2 + 1 constants c m , m = 0, • • • , n 2 and the q\ constants 
b( y £ = 1, • • • , qi . This raises the question of an analogous approach for 
q 2 > 1. We have carried out the corresponding analysis and will present 
the results in a separate paper. The results are useful, in particular, 
when gi is large, or even infinite. 

VII. PRIMARY QUEUE ONLY: CASE II 

We now begin the investigation of the second case, in which no 
overflow is permitted from the primary queue, so that a demand in the 
primary queue must wait for a server in the primary to become free. 
Hence, R\ 2 = in the diagram in Fig. 1, and this is the only aspect in 
which the second case differs from the first. In this section, we consider 
the case in which there is a primary queue only, with q\ > 1 waiting 
spaces. The steady-state probabilities p$ now satisfy the equations 

[ai(l - S iki Sjn 2 ) + a 2 (l - 8jn 2 ) + min(i, m) +j]pij 

« <*l(l - fi,- )[l - X*-rt,-l(l - 8jn 2 )]Pi-l.j 

+ (1 - 8jo)(aiXi-nj + a*)pij-i 

+ (1 - 5 lftl )min(i + 1, »i)p, +w + (1 - 8 Jai )U + l)Pu*i* (109) 
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for < i < ki, < j < /i 2 , where Ai - »i + gi. The normalization 
condition is 

But (109) for < i < ni - 1 and <y < n z leads, analogously to case 
I, to the representation 

n 2 

Po = S D m Si(p m , ai)8j(-p m , a 2 ), < i < m, <> < n 2 , (111) 

m=0 

where p m , m = 0, • • • , n 2 are the roots of (10), with p = 0. The 
coefficients D m are to be determined. Also, for n\ + 1 < i < Ai and 
0<y</i 2 - 1, 

(ai + a 2 + n\ + j)pij 

= (ai + a 2 )(l - 5,o)/>.j-i + /ii(l - hkjpi+ij + U + Dpv+i- (112) 

We define the quantities IT r7 , for r,j = 0, 1, • • • , as the solutions of 
the equations 

(ai + a 2 + ni +y)n ry 

- (a! + a 2 )(l - 8jo)U r j-i + «i(l - MU r -ij + (j + l)n r , y+ i, (113) 

which satisfy the initial conditions 

n r0 = 5 r0 , r - 0, 1, ... . (114) 

It follows from (163) that 

£ n ry y r = 8y(»i(l - y), d + a 2 ), (115) 

r-0 

and, in particular, 

n 0y = Sj(n u ai + a 2 ). (116) 

Also, (115) implies that 

n r , = 0, r>j. (117) 

It may be verified, from (112), (113), and (117), that 

*i 
Pij = I ProUr-u, ih+lSiS*i, <y < n 2 . (118) 

The quantities IIi-> may be calculated sequentially from the recurrence 
relations (113), with the help of (117) and the initial condition n o 
= 1. 

It remains to satisfy the boundary conditions at i = ri\ , < j < n 2 
and at j = n 2 , n\ + 1 < i < ki, as well as the normalization condition 
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(110). Hence, from (109) we must have 

[ai + a 2 (l - Sj n2 ) + ni +j]pn iJ = aip^-ij + (ai + a 2 )(l - 8jo)p ni j-i 

+ rupn^ij + (1 - 8j n2 )(j + l)p„„, + i (119) 
for < j < n 2 and 
[ai(l - fi.Aj) + /ii + n 2 ]A> 2 

= aiA-i,n 2 + (ai + a 2 )pi,n 2 -i + ni(l - SiA^p.+i,^ (120) 

for ni + 1 < i < £i . 

From (111), (118), and (119), it is found, after reduction, that 

£ D m [p m s ni (l + p m , ai)sj(-p m , a 2 ) + ais n ,(p m , ai)sj(-l - p m , a 2 )] 

m-0 

-m E Pron r _ ni - liy = (121) 

r=n, + l 

for <y < ri2 - 1. Also, with the help of (113), the boundary conditions 

(120) may be written in the form 

£ Pro[(il2 + l)n r - ni _i,„ 2 + i — (ai5n,+l,A, + 02)n r _« 1 -i,rt 2 ] 
r-71, + 1 

"2 

- ai £ DmSnSpm, ai)Sn 2 (-pm, d 2 ) = 0, (122) 
m— O 

and, for q\ > 2, 

S Pro[(n 2 + l)n r -,,n 2+ i - (ai5,Ai + a 2 )n r _,, n2 ] 

r— i 

-ai 2 p, n r - l+ i,„ 2 = 0, (123) 

r-i-l 
for Mi + 2 < j < ^i . 

The boundary condition corresponding to y" = n 2 in (119) may be 
reduced, with the help of (10), to that in (121) fory = n 2 . However, this 
last condition may be shown to be redundant, and it suffices to use 

(121) for <y < n 2 — 1. These conditions, together with those in (122) 
and (123), determine the constants D m , m = 0, •••, n 2 and the 
probabilities p r o, r = ni + 1, ■ ■ • , k\, only to within a multiplicative 
constant, which is determined from the normalization condition (110). 
But, from (111), with the help of (10) and (171), it follows that 

t Pij = DosAO, ai)s n2 {l, a 2 ), < i < n x . (124) 

7=0 

Also, from (109), if we sum ony from to n 2 , and on i from /to ki, we 
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obtain 



ni 2 Pa = aipt-i,^, Bi+ls/sAi. 



(125) 



7-0 



Hence, from (124) and (125), with the help of (111), (118), and (171), 
the normalization condition (110) implies that 

a\ " 2 
DoSnAl, ai)«u(l, a 2 ) + — 2 D m s ni (p m , ai)s„ 2 (-p m , a 2 ) 
n.\ m =o 



ai 



r— «i— 1 



+ - S ph) S n^-p^ono^l-i. (126) 

Once the constants Z/ m , m — 0, • • • , n 2 and the probabilities p r o, r 
— Tii + 1, • • ■ , ki have been determined, the remaining probabilities 
may be calculated from (111) and (118). Also, other steady-state 
quantities of interest may be calculated directly. In particular, the loss 
probabilities are 

L\ = Pk u n 2 = Pki.0U0.n2, (127) 

and, with the help of (171), 



2 

i=0 



"2 

2 

m=0 



"1 "2 

Ll= 2 P'.n 2 = 2 A»Sn,(l + Pm, ai)Sn 2 (-pm, 02) 



+ J /> r0 S -W (128) 

r=/i 1 + l <f-0 

The probability that a demand from the primary source is queued on 
arrival is 



"2 

2 

m=0 



re]— i .. z 

Ql = 2 P«'.»2 = 2 DmSn x (Pm, d^Sn^—pm, fife) 



Aj 1 — wj— 1 

+ x Pro 2 n^ 2 -p Al , n ,„ 2 , (129) 

r-fij+l ^=0 

which may be simplified in view of (126). 

The probability that a demand arriving from the primary source 
overflows (immediately) is 

A, n 2 -\ H2-1 / A, \ 

/i2= 2 2 pu= 2 [Pn 1 j+ 2 PiA- (130) 

»=n, y'-O y'=0 \ «=rt,+l / 

With the help of (125) and (171), it is found that 



/l2 = 2 DmSnApm, d\) 



d\ 

S*,-l(l — Pm, OI2) H S„ 2 (— p m> 02) 

Hi 

r— ti\— 1 



+ (— -1) 2 PrO 2 n^-— p^.ono.n,. (131) 
K R\ I r=n,+l <?=o n.\ 
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The mean departure rate from the primary queue to the primary 
servers is 

Rn-nip £ "ipij. (132) 

i=nj+l y=0 

It follows from (1) and (125) that Ru = AiQi, as is to be expected. 
The average number of demands in the primary queue is 

Vi- £ I (i-n 1 )p u . (133) 

i=n,+l y'=0 

Again with the help of (125), it is found that 






"2 

2 D m S ni (pm, Cli)Sn 2 (—pm, d 2 ) — (?1 + l)pk u oIlo,n 2 
m-0 

*1 r— nj— 1 

+ I Pro % (r+l-Z-nOn** 



(134) 



The average number of demands in service in the primary is 

Xi- J 2 minfe in)**. (135) 

i=0 7=0 

With the help of (110), (124) and (172), we obtain 

Xi - »i - D s„,-i(2, ai)s„ 2 (l, 02). (136) 

The average waiting time of the demands which are queued in the 
primary is W\, as given by (87). Also, since i?i2 = and Xi = aiu, the 
first relationship in (88) becomes 

ai(l-L,-/ 12 )=A-,. (137) 

In view of (127), (131), and (136), this relationship provides a useful 
numerical check. Finally, we comment that the only U. r j which occur 
in (127) to (129), (131), and (134) are U^, for <f= 0, • ■ • , q x - 1. 

VIM. PRIMARY AND SECONDARY QUEUES: CASE II 

We now consider the second case when there is a secondary queue, 
with q 2 > 1 waiting spaces, as well as a primary queue, with <?i > 1 
waiting spaces. The steady-state probabilities satisfy the equations 

[a,(l - SikJtj-nJ + a 2 (l - 8jk.) + min(i, m) + mm(j, n 2 )]p,j 

= d(l - 5,o)(l - Xi-n,-\X" 2 -\-j)Pi-\J 

+ (1 - 8j< ) )(a i x>-n ] X">-j + a ^P'J-i + (1 - 5 «0 min (* + !. »i)p*fU 
+ (1 - 8/*a) min(y + 1, n 2 )pu+h (138) 

for < i £ ki, </' < & 2 , where £i and £ 2 are as defined in (59). The 
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normalization condition is 

*l *2 

S ZP.;=1- (139) 

1=0 y=0 

Now (138) for < i < /ii — 1 and < j < k 2 leads, analogously to 
case I, to the representation 

k 2 

£ C m Si(p m ,ai)sj(-p m ,a2)6 ni (pm)<fyn. 2 (pm), 0</<n 2 , 

, m-0 

Pu-\ k < 140 > 

2 C m Si(p m , ai)s n ,(-p m , a 2 )0 ni (p m )(f>Apm), n 2 <j<k 2 , 

m-0 

for S i S Hi, where p m , m = 0, • • • , k 2 are the roots of (30), with po 
= 0, and the constants C m are to be determined. Similarly, (138) for 
m + 1 < i < Ai and n 2 + 1 <,/ < k 2 , analogously to case I, leads to the 
representation 

* 2 

Pij= £ C m S ni (Pm, Cli)Sn 2 (—pm, a 2 )di{p m )<t>j(p m ) 
m-0 

Q\ 

+ £ B £ 6i(v,)4>j(vt), ni < i < A,, n 2 <y < £ 2 , (141) 
e-\ 

where v,, <f = 1, • • • , q\, are the roots of On x (v) = 0, and the constants B, 
are to be determined. As before, the representations in (140) and (141) 
agree for i = tli. 

Finally, (138) for /ii + 1 < i < k\ and <y < >i 2 — 1 leads, as in the 
previous section, to the representation 

A, 

Pii = 2 ProUr-u, Hi + 1 < i < M, < y < « 2 . (142) 

In order that the representations in (141) and (142) agree for j = n 2 , we 
must have 

A, *2 

X ProTIr-j^ = J C m S ni (pm, ai)S^(~pm, a2)^(pm)<f'n 2 (P'") 



+ I BAM^v,), m + 1 < i < *i. (143) 

Since, from (116) and (167),n ,„ 2 = s„ 2 (ni, ai + a 2 ) > 0, these equations 
may be used to solve successively forpA„o, • • • ,p Rl +i,o in terms of C m , 
ra = 0, ••- ,k 2 and B f , <f = 1, • • • , q\. 

It remains to satisfy the boundary conditions at i — Hi, £ / < Aa 
and at y = /i 2 , ni + 1 < i < &i, as well as the normalization condition 
(139). Hence, from (138), in addition to (119) for <> < n 2 - 1 and 
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(71) for n 2 + 1 £7 £ k 2 , we must have 
(ai + a 2 + ni + n 2 )p ny n 2 

= aiPn,-i,n2 + («i + a 2 )pn l ji 2 -i + niPm+i^ + n&n^+i, (144) 
and 
[ai(l - 6a,) + a 2 + ni + n 2 ]p,,„ 2 

= a1p.-1.74 + (ai + a 2 )pi^-i + /ii(l — 8 i k)pi+\^ 1 + n-iPi^+i, (145) 

for ni + 1 < i < Ai. 

From (119) for <y < n 2 - 1, (140) and (142), it is found, analogously 
to (121), that 

2 C m [p m s„,(l + p m , ai)s 7 (-p m , a 2 ) 

m-0 

+ ais„,(pm, ai)s>(-l -p m , a 2 )]^,(pm)^(pm) 

-ni 2 Pr n r - ni _ w = 0, (146) 

r=n,+l 

for < 7 < n 2 - 1. Also, from (71), (140), and (141), it is found, 
analogously to (74), that 

*2 

2 CmPmSn 2 (-pm, fl 2 )|> ni (l + Pm, ai)0„,(p m ) 
m-0 

+ S ni (pm, ai)tiq } -l(p m )]<i>j(pm) 

= m 2 BA 1+ i(^M, (147) 
<r=i 

for /i 2 + 1 <y < k 2 . From (141), (142), and (145), with the help of (177) 
and (180), it follows that 

a 2 £ CmS ni {pm,ai)Sn 2 {-p m , a2)d i (p m )^n. 1 -\{pm) 
m=0 

+ a 2 2 BAM4>n 2 -iM - (a, + a 2 ) £ ProIIr-4-,-1 - 0, (148) 

^=1 r-i 

for ni + 1 < i < fa. 

The boundary condition (144), which may be reduced with the help 
of (30) to that in (146) for,/ = n 2 , can be shown to be redundant. The 
conditions in (143) and (146) to (148) determine the constants C m , m 
= 0, • • • , k 2 and B,, /= 1, • • ■ , Oi, and the probabilities p,o, r = /ii + 1, 
• • • , fa, only to within a multiplicative constant, which is determined 
from the normalization condition (139). But, from (140), with the help 
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of (27), (30), and (171), it follows that 
2 Pu = CosdO, ai)tf Bl (0)[8^(l, a 2 )^ 92 (0) - s^-id, a 2 )% 2 -i(0)], (149) 

for < i < Tii. Also, from (138), if we sum on/ from to &2, and on i 
from «f to Ai, we obtain 

»i 2 P* - «i 2 P«j» ni + 1 < /< ki. (150) 

y'-o y'=i2 

Hence, from (149) and (150), with the help of (27), (66), (141), and 
(171), the normalization condition (139) implies that 

CoSn.U, ai )e ni (0)[sn 2 (l, a 2 )% 2 (0) - s^-iU, a 2 )% 2 -i(0)] 

ai * 2 

+ — 2 GnSn,(Pm, ai)s^(—p M , a 2 )[tt qi (p m ) ~ l]*q 2 (Pm) 
fll m _o 

+ - 2 BJ&S.P*) - i]%M) = i. (15D 

Once the constants C m , m = 0, • • • , k 2 and 5 A ^=1, • • • , qi, and the 
probabilities p,o, r = n\ + 1, • • ■ , k\ have been determined, the 
remaining probabilities may be calculated from (140) to (142). 

IX. SOME STEADY-STATE QUANTITIES: CASE II 

Other steady-state quantities of interest are depicted in the diagram 
of Fig. 1. The loss probabilities L\ and L 2 are given by (78), and the 
queuing probabilities Qi and Q 2 are given by (79). Also, the mean 
departure rate from the secondary queue, R 22 is given by (82). The 
average number of demands in the secondary queue, V 2 , is as given by 
(84), but the average number of demands in the primary queue is 

V, = 2 2 (* " ndPv- (152) 

i=n,+l >=0 

The average waiting times of the demands which are queued in the 
primary or in the secondary, W\ and W 2 , are given by (87). The 
(immediate) overflow probability I i2 is given by (130). The mean 
departure rate from the primary queue to the primary servers is 

ft] k 2 
Uii-nifi 2 2P«y> (153) 

and the average number of demands in service in the primary is 

X, = 2 2 min(i, n,)p iy , (154) 

,=o y=o 

The above steady-state quantities may be expressed in terms of the 
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constants C,„ and B,, with the help of the representations in (140) and 
(141) and the relationship (150). Because of (150), we do not have to 
use the representation in (142), and hence we are able to eliminate the 
quantities 11,-/. From (78), analogously to case I, and corresponding to 
(89) and (90), 

*2 </l 

L, = £ C m S ni (pm, adSn^-pm, a 2 )* q2 {f>m) + £ B& q2 {Vf) (155) 
m=0 ('-l 

and 

*2 

L 2 = X C m S n . 2 (-pm, a 2 )[Sn,(l + pm, di)6„ x (p m ) 
m-0 

+ S ni {Pm, aJQ^-dprn)] + £ IWVlOv). (156) 

Similarly, from (79) we obtain, corresponding to (91), 

*2 
Ql = 2 CmSn^Pm, d\)Sn^-pm, a 2 )[ft 9l (p m ) ~ l]¥<7 2 (pm) 
m-0 

+ 2 ii^ - rPfefc), (157) 

which may be simplified in view of (151). 

From (130), with the help of (27), (66), (150), and (171) it is found 
that 

*2 
/l2 = £ CmSn^pm, a])S„ 2 -l(l - p m , a^tfn^Pm^n^pm) 
m-0 

*2 

+ X CmS n ,(Pm, d^Sn^—Pm, d 2 ) 

m-0 

- [a,,(p m ) - i] - Qn-ifoj 1 *,>«) 

ni J 

^=i l.Wi J (158) 

From (1), (150), and (153), it follows directly that R n = XiQu as is to 
be expected. From (154), with the help of (139), (149), and (172), it 
follows that 

Xi - ni - Cos„,-,(2, a,)0„,(O)[S",(l, d 2 )* g . 2 (0) 

- s„ 2 -id, a 2 )¥„ 2 -i(0)]. (159) 

In view of (155), (158), and (159), the relationship in (137) provides a 
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useful numerical check. From (152), with the help of (27), (51), (66), 
(150), (182), and (183), it is found that 

Vi = C o s„,(0, a l )s ri2 (0, a 2 )* 92 (0)A 9l I— 

+ — Yi C' m s„ 1 (p m , ai)s„ 2 (-p m , a 2 )^ 92 (pm) 

fl\ m = \ 

.{«,>«) - (qi + 1) - — [1 - Or,(P*) + -fi, r i(pJ] j 

l_ Pm ai J 

+ - I B,% 2 M K,(zv) ~ (<7> + 1) 

- - [1 - Sl qi M + - Sl^b*,)]). (160) 
Vt ai J 

Now (60) and (138) are identical for < i < k\, n 2 + 1 < y < & 2 . 
Hence, if we define r y as in (98), and proceed as in case I, we deduce 
that (50) and (52) hold, and R 22 = A2Q2, as is to be expected. 

Although we assumed here that q 2 > 1, we remark that the expres- 
sions in (155) to (160) remain valid for q 2 = 0. This may be verified 
from (127) to (129), (131), (134), and (136), by substituting D m = 
C m ni (pm) and using (143), noting that^ 2 (p) = 1, and k 2 = n 2 for q 2 = 
0. We point out that 

t P* r "i n,„ 2 = i E^n r _,„ 2 , (i6i) 

r=n, + l /=0 i=/i, + 1 r=i 

and 

£ Pk) I (r + 1 - €- n,)U,n 2 



= 2 (i + 1 - m) 2 PK.n r _ I>2 . (162) 

Similarly, the normalization condition (151) holds for g 2 = 0, as may 
be verified from (126). Also, the boundary condition (146) reduces to 
(121) for q 2 = 0. Finally, it may be verified, with the help of (113) and 
(143), that the boundary conditions (122) and (123) are equivalent to 
(148) for q 2 = 0. 
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APPENDIX A 

We define s,(A, a) by the recurrence relation 

(a + i + A)s,(A, a) = a(l - 5, )s,_ 1 (A, a) + (i + l)s,+i(A, a); 

so(A, a) = 1, (163) 

for i = 0, 1, • • • . Thus s„(A, a) is a polynomial of degree n in both A and 
a, and it may be related to a Poisson-Charlier polynomial. 7,8 However, 
we give here the properties of s„(A, a) which we will need. It follows 
from (163) that the generating function 

g(u) = £ Si(X, a)u l (164) 

1=0 

satisfies the equation 

(1 - u) ^- = [a(l - u) + X]g; g(0) = l. (165) 

Hence, 

£ s,(A, a)u l = e""(l - u)-\ (166) 

1-0 

This leads to the explicit formula 

8i(X, a) = £ -|^t, (167) 

^o r\ (i-r)\ 

where 

(A) = i, (A) r = A(A + 1) • • • (A + r - 1), r = 1, 2, . • • . (168) 

If we differentiate (166) with respect to u, we obtain 

(i + l)s, + ,(A, a) = as,(X, a) + X s,(A + 1, a). (169) 

Also, if we write (1 - u)~ K = (1 - u)(l - u)~ x '\ it follows from (166) 
that 

s,(A, a) = s,(X + 1, a) - (1 - 5,o)s,-,(A + 1, a). (170) 

From (170), we obtain 

£ s,(A, a) = s„(X + 1, a). (171) 

From (169) and (171), we deduce that 

£ (n - r').s,(A, a) = (1 - 5„o)s„-,(A + 2, a). (172) 

/=0 
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We now turn our attention to the Chebyshev polynomials of the 
second kind, 9 UAx). They may be defined by the recurrence relation 

2x UAx) - U,+i(x) + Ut-i(x); £/_i(*) ■ 0, U (x) s 1, (173) 

for £= 0, 1, • • • . From (26) and (173), it follows that 

(a 2 + n 2 - p)*Ap) - 02 ¥*i(p) + i 2 ^i(p); 

*_,(p) - 0, * (p) - 1. (174) 

Since 9 



Ui 



it also follows that 



iH 



<r/2 



= 2 *~, 

r=0 



MM... -I*) Ui* + * -2(5,. 

• a 2/ \2Va 2 /l2/ r-0 \«2 y 



(175) 



(176) 



From (27) and (174), we deduce that 

[a 2 (l - 5/Aj) + »2 - p]4>/(p) = at&j-Ap) + n 2 (l - Sjkjto+iip), (177) 

for/ < ^2- 

Similarly, from (65) and (173), it follows that 

(ai + «i + p)SIAv) = ai&lfin(v) + HiQ^-i(r); 

Q-i(ir) ■ 0, Qo(f) - 1. (178) 

Also, from (175), we obtain 



m 



.«!/ \2VoiW,/ r=o\a.' 



(179) 



From (66) and (178), we deduce that 

[ ai (l - V + ii + v\Bi(v) m mOi-Ap) + »i(l - 8 ikl )e i+1 (v), (180) 
for i < ki. Also, from (179), 



fc(O) = 



Hence, since &i = /ii + oi, 



J (i , -B,)ft(0)-A, l -), 

i-n,+l \Oi 



ft] 



(181) 



(182) 



where A 9 (£) is as defined in (51). If we multiply (180) by (i — n{) and 
sum on i, we obtain 
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v X (i - rn)Oi(v) = a, £ Bi(v) -m X 0,iv) 

i=n, + l «=»i i-n,+l 

= aiEOg.W - 1] - inO^-idO, (183) 
from (66). We also note, from (178), that 

aift»» - n,6 ni+i (p) = iQn-iCF). (184) 



APPENDIX B 

We first consider the roots of s„(l - p, a) = 0, where a > 0. We 
relate this to an eigenvalue problem for a real symmetric matrix, so 
that we may apply a result of Slepian and Landau 12 to prove that the 
roots are distinct. Fory = 0, 1, • • • , we let 



J- 
yj= \/-3«>(i -p»«)- 



Then, from (163) it follows that 

v5>y-i + -=(p-l-a-j)yj+ y/j+lyj+i = 0. 
va 

If s„(l - p, a) = 0, then y„ = and 



(yo, • • • , y n -i) T = 0, 
where T denotes transpose, I„ is the n X n unit matrix, and 



A„ ■= (1 + a - p)In 

va 



A, = [0], A„ = 



1 
va 



vn- 1 - 



(n-1) 

Va" 



n>2. 



(185) 



(186) 



(187) 



(188) 



Hence the equation s„(l - p, a) = is equivalent to 



det 



A„ - -= (1 + a - p)\ n 
va 



= 0. 



(189) 



We remark that Kaufman, 6 who has analyzed the same overflow 
problem using matrix techniques, arrived at the same eigenvalue 
problem by block diagonalization. 

Now A„ is a real symmetric matrix, so that its eigenvalues are real. 
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Moreover, A„ may be obtained from A„+i by deleting the last row and 
column. Hence we may apply the results of Slepian and Landau. 12 In 
order to establish that the eigenvalues of A„ separate those of A„+i for 
n = 1, 2, • • • , we have to verify that certain quantities (see (6) of Ref. 
12) are different from zero. Specifically, let a be the n-dimensional 
column vector obtained from the last column of A„+i by omitting the 
diagonal element, so that 

r u),/i = i, 

I (0, ...,0, y/n) T , n>2. (190) 

Hence 

(a, (yo, ••-, y n -i) T ) = Jny n -i * 0, (191) 

since y n -i ¥= for any eigenvector corresponding to y n = 0. This 
establishes the desired result, and we may conclude that the zeros of 
s„(l — p, a) separate those of s„+i(l — p, a), for n = 1, 2, • • • . In 
particular, the zeros of s„(l — p, a) are real and distinct. Since, from 
(167), s„(X, a) > for X > 0, it follows that the zeros of s„(l - p, a) are 
larger than 1, and hence positive. 
We next consider the zeros of 

PnM = s n (l - p, a)\p q (p) - s„-i(l - p, a)\p q -i(p), (192) 

where 

We let 




z J=-\ /^_ L ^_.\ (194) 

j = n - 1, ■••,n + q, 

where yj is as in (185). We note that z n + g = 0, since U-i(x) = 0. Also, 
from (185), (192), and (193), the consistency of the definitions of z n -\, 
and of z n , in (194) implies that P n , q (p) = 0. From (186), it follows that 

yfj Zj-i +—={p-l-a- j)zj + V/+ 1 Zj+i = 0, 

Vfl 

j = o, . . . , n - 1, (195) 
and from (173) we obtain 

yfn Zj-i + —= (p — a — n)zj + Vn zj+i = 0, 

j = n, • • • , n + q - 1. (196) 
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But P n ,o(p) = s n (l - p, a) and we have shown that its zeros are real and 
distinct. Proceeding as before, it follows that the zeros of Pn.gip) 
separate those of P„. v+ i(p), for q = 0, 1, • • • . In particular, the zeros of 
Pn, q (p) are real and distinct. 

We now consider p < 0, and will show that P n , g (p) i* 0, for n = 1, 2, 
. • . and q - 1, 2, • ■ • . Since, from (170), 

s„(l - p,a) = s n (-p, a) + s„_i(l - p, a), (197) 

and, from (167), s„(A, a) > for X > 0, it follows that s n (l - p, a) > 
s n -i(l - p, a) > for p < 0. Hence 
s„(l - p, a) 



s„-i(l - p, a) 
Also, for p < 0, 



>1, p<0. (198) 



a + n-p 1/ , 1\ „_ //* /« 



-iU + - , ^maxU ,W- (199) 



But, from (175), 



9 

„2r-q 



•nU^iid^Hr, + (1/ti))] _£ <1 (2Q0) 

t/,[(l/2)fo + (1/tj))] £ ^ 

It follows from (193), (199), and (200), that 

^-i(p) .1 /a 



<- A /-<l, p<0. (201) 

hip) I Vn 

From (192), (198), and (201), we deduce that P n , q (p) 5^ for p < 0. 
Thus the zeros of P n , q (p) are positive. 

Finally we consider the zeros of 0„, (v). For <f = 0, 1, 2, • • , we define 

wM , V i«±*+1\ -J«Vn (* + £ii) . (202) 

Then, from (59), (65), and (66), n ,(p) = is equivalent tou^OO = 0. 
From (173) we deduce that 

(/ii + p)u>o = "Ja\n\ w\, (203) 

and 

Vai/ii w/-i - (ai + Aii + r)«fc + Vaiiii it?*.] = 0, (204) 

for <?= 1, 2, ■ • ■ . Proceeding as before, it follows that the zeros of w q (v) 
separate those of w q +\(v), for q = 1, 2, • ■ • . In particular, the zeros of 
w q (p) are real and distinct. But, from (193), (201), and (202), it follows 
that w q (v) ¥" for v > 0. Hence, the zeros of w g (v) are negative. 
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